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Abstract 

We extend the definition of the Costas property to functions in the continuum, namely on intervals of the reals 
or the rationals, and argue that such functions can be used in the same applications as discrete Costas arrays. 
We construct Costas bijections in the real continuum within the class of piecewise continuously differentiable 
functions, but our attempts to construct a fractal-like Costas bijection there are successful only under slight 
but necessary deviations from the usual arithmetic laws. Furthermore, we are able, contingent on the validity 
of Artin's conjecture, to set up a limiting process according to which sequences of Welch Costas arrays converge 
to smooth Costas bijections over the reals. The situation over the rationals is different: there, we propose an 
algorithm of great generality and flexibility for the construction of a Costas fractal bijection. Its success, though, 
relies heavily on the enumerability of the rationals, and therefore it cannot be generalized over the reals in an 
obvious way. 

1 Introduction 

Costas arrays [3] have been an active topic of research for more than 40 years now; however, after 1984, when 2 
algebraic construction methods for Costas arrays were published (the Welch and the Golomb method [7]), still 
the only ones available today, there has been effectively no progress at all in the construction of new Costas 
arrays, with the obvious exception of brute force searches. Recent research on Costas arrays tends to focus on 
the discovery of new properties [HO [10], hoping that they will either furnish some lead for a new construction 
method, or prove that such a method does not exist, and thus overcome the current virtual stalemate in the 
core problems of the field. 

In line with this effort, it is likely that research on Costas arrays would benefit by the extension of the 
definition of the Costas property in the continuum, for 2 reasons: on the one hand, this might open the door to 
assistance from the entire arsenal of analysis, as was the case with the successful generalization of the factorial 
in terms of the Gamma function; on the other hand, the recent advances in the subject of the Instantaneous 
Frequency of a signal [H] make it possible to design signals with continuously varying frequencies instead of 
piecewise constant frequencies, such as the usual discrete Costas arrays model, and there might be benefits in 
doing so. And besides, such objects certainly have an intrinsic pure mathematical merit for study. 

In this work, we propose a suitable extension of the definition of the Costas property in the continuum (which 
we take here to mean the real and rational numbers), and we explain how the existing discrete Costas permuta- 
tions can be used to generate continuum Costas permutations. Note that, in accordance with common practice 
in recent literature, we will be using the terms "Costas permutation" and "Costas array" interchangeably. 

*The author holds a Diploma in Electrical and Computer Engineering from NTUA, Athens, Greece, and a Ph.D. in Applied and 
Computational Mathematics from Princeton University, NJ, USA. He was a scholar of the Lilian Boudouris Foundation, 
t www.shannoninstitute.ie 
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2 Basics 

We reproduce below the definition of a Costas function/permutation [3]: 

Definition 1. Let [n] := {0, . . . , n — 1}, n £ N and consider a bijection / : [n] — > [n]; / is a Costas permutation 
iff the multiset {(i — j, f(i) — f(j)) : < j < i < n} is actually a set, namely all of its elements are distinct. 

These permutations are extremely useful because they give rise to binary signals with an optimal autocor- 
relation pattern: 

Definition 2. Let / : [n] -> [n], n 6 N* be a Costas permutation, and let F : 1? — » [2], the corresponding 
binary signal of /, satisfy F(i, f(i)) = 1, i £ [n], and F — everywhere else. The autocorrelation of / is: 

A F (u,v)= F(u + i,v +j)F(i, j), (u,v)£Z 2 

The following result is just a restatement of the Costas property: 

Theorem 1. Let f : [n] — > [n], n £ N*, be a permutation, and let F be its corresponding binary signal; then, 
< Af(u, v) < 2, Vu, v £ Z 2 — {0, 0} ij(f / /las the Costas property. 

We have already mentioned the Welch construction method for Costas arrays. As we will refer to it several 
times below, we offer its definition for the sake of completeness: 

Theorem 2 (Welch construction Wi (p,g,c)). Let p be a prime, let g be a primitive root of the finite field ¥(p) 
of p elements, and let c £ [p— 1] be a constant; then, the function f : [p — 1] + 1 — > [p— 1] + 1 where f(i) — </ 1+c 
mod p, i £ [p — 1] is a bijection with the Costas property. 

3 Costas bijections in the real continuum 

From now on, until Section [7] we will be using the term "continuum" in the sense of "real continuum", unless 
explicitly stated otherwise. 

3.1 Definitions and simple results 

In our extension of Definition[T]in the continuum we will replace [n] by [0, 1], but otherwise the definition remains 
the same: 

Definition 3. Consider a bijection / : [0, 1] — > [0, 1]; / is a Costas permutation iff the multiset {(x — y, f(x) — 
f(y)) '■ < y < x < 1} is actually a set, namely all of its elements are distinct. 

Remark 1. The choice of the interval [0, f] is by no means restrictive: it can be seen immediately that for any 
pair a, b £ R, a < b there exists a linear monotonic mapping h mapping [0,1] bijectively on [a, b], specifically 
h(x) = a + x(b — a), < x < 1, and / has the Costas property on [a, b] iff /i -1 o / o h has the Costas property 
on [0, 1]. 

Yet again, we can give an alternative but equivalent definition of the Costas property in terms of autocorre- 
lation: 

Definition 4. Consider a bijection / : [0, 1] — > [0, f], and let F : M 2 — > {0, 1, oo} be its corresponding quasi- 
binary signal (that is, binary whenever finite), so that F(x,f(x)) = 1, x £ [0, 1], and F = otherwise. The 
autocorrelation of / is: 

Af(u,v)= f f 5(F(x + u,y + v) - F(x,y))dxdy, (u,w)£R 2 
Jo Jo 

Remark 2. Notice that this autocorrelation, just like its discrete counterpart in Definition [2] takes integer 
values whenever finite, as it counts the number of zeros in the argument of the Dirac ^-function. 

Once more, then, the following result is just a restatement of the Costas property: 

Theorem 3. Consider a bijection f : [0, 1] — ► [0, 1], and let F be its corresponding quasi-binary signal; then, f 
has the Costas property iff < Af(u,v) < 2, V(u,v) £ R 2 - {0,0}. 
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3.2 Applications 

Continuum Costas bijectionfl can find applications in the same situations their discrete counterparts do [3]- 
For example, consider a RADAR system whose operation relies on a usual Costas waveform. In practical terms, 
this means that the waveform it transmits is of the form: 

w(t) = ylcos j 2iv j - - — /lri T fc+i T \ (t) j t J , s a Costas permutation of order n, t £ [0,T) 



/ s(k) \ 

which is a different way to express that, for n G N*, w(t) = A cos ( 27r— — ft I , t €E 

V n J 

Alternatively, we could have used a continuum Costas permutation s on [0, 1] . Let us consider the waveform: 



t T ,tt± T ), ke[n] 
n n 



w(t) = A cos ^2vr/ ^ s(w)dti + 2tt/ ^ , * € [0, T) 
Bedrosian's theorem [HE] on instantaneous frequency asserts that the instantaneous frequency of w is 

(lirf J* s(u)du + 2irf tj = s(t)f + f , 
as long as w(0) = 0; this condition can be satisfied, at least approximately, through an appropriate choice of /q. 



3.3 Link between continuum and discrete Costas permutations 

How do the 2 definitions compare? The expression for the discrete waveform is clearly a special case of the 

s(k) + 1 

continuum expression, and this can be seen if we write S(t) = /Ifiy k+i T \(t), where s is a Costas 

k = 

array of order n and S is a continuum permutation (but obviously not Costas). The verification of the Costas 
property through the autocorrelation in the discrete case is also a subprocess of the verification in the continuum 
case: we just need to take care that horizontal and vertical displacements of the copies of the functions in the 

T f 

autocorrelation formula are integral multiples of — and — , respectively. 

n n 

As S is a piecewise constant function, one might be tempted to attempt to formulate a definition for (at least 
a class of) continuum Costas permutations in terms of Costas arrays, as limits of sequences of Costas arrays, 
just like measurable functions are approximated by sequences of piecewise constant functions: a Costas array 
s n of order n can be mapped on a piecewise constant function S n , just as we did above, and, letting n — > oo, we 
can hopefully obtain a continuum Costas permutation S. This limit would probably be highly discontinuous, of 
a fractal nature perhaps, as Costas arrays are highly erratic and patternless. 

The problem with the plan of action suggested above is that we seem to have no good understanding yet of 
sequences of Costas arrays across different orders that follow a clear pattern, so that we can successfully describe 
how the limit of such a sequence would look like; a notable exception is the example we give below in Section [6] 
Nevertheless, the idea of seeking continuum Costas permutations among fractals seems, in principle, promising 
in itself and worthwhile investigating. But first, let us focus on the case of smooth functions. 



4 Construction of smooth continuum Costas permutations 

The whole idea of the existence of smooth functions with the Costas property may sound outright irrational 
at first, and any investigation futile: after all, there can hardly be any object more irregular and discontinuous 
that Costas arrays. Nonetheless, the continuum is dense in itself, while finite discrete sets are not, and this 
makes a big difference, as we are about to see: for example, the function fix) = x 2 has no chance of being a 
permutation on any discrete set other than [2] U {oo}, while it is a permutation on both [0, 1] and [1, +oo], as 
it effectively makes some areas of the intervals "denser" and some "sparser" (consider, for instance, that the 
images under / of all points in [0, y2/2] get "crammed" in the smaller interval [0,0.5]). In the continuum we 
can create Costas permutations by causing "elastic deformations", by "changing the density" of points in an 
interval, whereas such techniques are inapplicable on discrete sets. 

1 We have to resort to the use of the uncommon word "continuum" in the role of an adjective here instead of the perhaps more 
appealing intuitively "continuous": the term "continuum function" accurately describes a function defined on an interval, or on 
something non-finite and dense at any rate, whereas the term "continuous function" has an already established different meaning in 
mathematics. 
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Let us begin by seeking functions with the Costas property that are reasonably smooth; for example, let us 
confine ourselves to special categories of almost everywhere differentiable bijections. 

Definition 5. Let / : [0, f ] — ► [0, 1] be a bijection; 

• / will be piecewise continuously differentiable iff there exists n £ N* U {oo} and a sequence of intervals 

o 

{Ii}2=i so that, for each i = 1, . . . , n, / is continuously differentiable in Ii (n = oo is used as a convention 
to denote a countable infinity of intervals); 

o 

• if, in addition to being piecewise continuously differentiable, for each i = 1, . . . ,n f \Ii is strictly monotonic, 
/ will be called piecewise strictly monotonic piecewise continuously differentiable; 

• if, in addition to being piecewise continuously differentiable, / satisfies the property that, for all sequences 

o 

of points such that Xi £ Ii, i = 1, . it is true that the sequences {/'(iCi)}™-! are either all 

strictly increasing or all strictly decreasing, / will be called overall strictly monotonic piecewise continuously 
differentiable; 

• / may combine all 3 features above, in which case it will be called overall and piecewise strictly monotonic 
piecewise continuously differentiable. 

Theorem 4. Let f : [0, f] — > [0, f] be an overall and piecewise strictly monotonic continuously differentiable 
bijection. Then, f has the Costas property on [0, 1]. 

Proof. Let us choose 4 points in [0, f], say x, y, x + d and y + d so that y < x and d > 0; these may actually be 
3 points if x = y + d. We need to show that 

f{x) - f{x + d) = fiy) -fiy + d)^d = 

Exactly one of the 2 pairs of intervals [x, y], [x + d, y + d] or [x, x + d], [y, y + d] consists of intervals with disjoint 
interiors. Without loss of generality, assume it is the second pair, then the Newton-Leibnitz Theorem implies 
that ^ ^ 

fix + d) - fix) = f f'iu)du, fiy + d)- fiy) = f f'(u)du 

J x J y 

Now, if / is overall and piecewise strictly monotonic continuously differentiable, it is always the case that either 
Vu, v : u £ ix,x + d), v £ (y,y + d)f\u) < f'iv) or Vu, v : u £ ix,x + d), v £ iy.y + d)f'iu) > /'(«), so that 
fix) - fix + d)^ fiy) - fiy + d) unless d = 0. □ 

Theorem 5. Let f : [0, f] — > [0, f] be a piecewise continuously differentiable bijection; if /' is not injective, f 
does not have the Costas property. 

Proof. We distinguish the following cases: 

• /' is constant on an interval, say /' = c € K , or, equivalently, / is linear on that interval: it follows there 

■ < ■ , , , • , , , „ , fix + d)- fix) fiy + d)- fiy) 

exist 4 points x, y, x + a and y + a with y < x and d > so that — = — = c, 

d d 

hence the Costas property is violated. 

• Assume that /' is never constant on an interval. Then, either there exist ii, i% so that \ f'ih 1 ) fl /' (7; 2 ) | > 0, 
namely it fails to be overall strictly monotonic, or there exists an i for which f'\Ii is not monotonic. In 
either case, there exist 2 points X\,X2 £ (0, f), so that xi < xi and f'ixi) = /'(a^). We distinguish 2 
subcases: 

— Neither of the points is an inflection point, that is both points lie in regions of the domain where / is 

either convex or concave; these regions are necessarily different, or the derivative could not possibly 

be equal at these points. This implies that there exist real numbers ei,62 > so that, if 2 parallels 

are drawn to the tangent at each of the points x\ and X2, at the side of the tangents where the 

function graph lies, and whose distances from the tangents are less than e± and ti, respectively, they 

each intersect the function graph at 2 points, say in < x\i. and £21 < 2:22- Clearly both x\\ — X12 

and X21 — X22 go to as the parallels move closer to the tangents, whence /(x'n) — /(112) and 

/(*2i) — fix 22) also go to 0; moreover, if ei and £2 are sufficiently small, (£11,2:12) n (0:21, X22) = 0, 

and each of (xn,a;i2), ix2i,X22) falls entirely within one of the intervals {Ii}, i = l,ldots,n. Hence, 

i ■ t 11 1 i-u + _ /(^la) fix2i)-fix 2 2) , 

we can choose a pair of parallels so that — = — and in —2:12 = X21 — X22- 

xu — X12 X21 — X22 

This violates the Costas property. 
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— At least one of the points is an inflection point, say xi, so there is a 5 so that x £ (xi — S,xi+S) — {xi} => 
f'(x) < ,f'{xi) and (ii — 8, x\ + 8) falls within one of the intervals {Ii}, i = 1, Idots, n, say Ik- As /' is 
continuous within Ik , and is not constant in any interval, there exist u\ £ (x\ — 5, xi) , u<x £ [x\ , x\ + 5 
so that neither is an inflection point and that /'(ui) = f'{u2). We are now back to the case above. 

□ 

Note that the derivative of a continuously differentiable bijection must keep the same sign throughout its 
domain, or else the bijection would have an extremum and would not be a bijection. Further, in the case 
of a continuously differentiable bijection, overall and piecewise strict monotonicity are identical, hence strict 
monotonicity implies injectivity. Therefore, in this special case, the following holds: 

Corollary 1. 

• Let / : [0, 1] — ► [0, 1] be a bijection continuously differentiable in (0, 1); then, / has the Costas property iff 
/' is strictly monotonic. 

• A continuously differentiable bijection on / : [0, 1] — * [0, 1] with the Costas property must be strictly 
monotonic. 

Remark 3. The issue of the continuity of the derivative of a function is quite esoteric. When a function is 
differentiable in an open interval, its derivative is not necessarily continuous. However, it is "almost" continuous, 
in the sense that, for any value between 2 values the derivative actually assumes at 2 points, there is a point 
between the 2 aforementioned points where the derivative assumes the chosen value. This property is known as 
Darboux continuity in the literature [2]. Working with piecewise continuously differentiable functions, we "float 
over" this technical point. 

Let us now see some examples of continuously differentiable bijections with the Costas property as well as 
some rules to produce new ones from known ones: 

Corollary 2. The following continuously differentiable bijections / : [0, 1] — > [0, 1] have the Costas property on 
[0,1]: 

• f(x) = x a , a £ R+, a ^ 0,1; 

• /(x) = ^— p a£R;-{l}; 

• f(x) = sin {^x 

Further, if /, g : [0, 1] — > [0, 1] are continuously differentiable bijections and have the Costas property on [0, 1], 
the following functions also do: 

• i-/; 

• af + bg, a, b £ R+, a + b = 1, if /, g are both strictly increasing or all strictly decreasing, and so are /',<?'; 

• fog, if f',g' are strictly monotonic of the same type and g is strictly increasing; 

• /ffj if /iSi/'jff' are au strictly increasing or all strictly decreasing. 

/ a x - 1 \ ' a x 

Proof. Observe that ( ) = ln(a) is strictly increasing for a > 1 and strictly decreasing for a < 1, 

\ a — 1 J a— 1 

(x a )' = ai° _1 > is strictly increasing when a > 1 and strictly decreasing when < a < 1, and ^sin — 

7T / TV \ 

— cos y~2 X j is strictly decreasing. Moreover, all of these functions are bijections, hence, they have the Costas 
property. 
Further, 

• (1 — /)' = — /' is strictly monotonic iff /' is, although of the opposite type, and 1 — / is a bijection on 
[0, 1], so it also has the Costas property. 

• (af + bg)' — af' + bg' is strictly monotonic if f' ,g' are both strictly monotonic of the same type, and 
af + bg is strictly monotonic too, hence a bijection, if /, g are both strictly monotonic of the same type. 

• / o g is clearly a bijection if both / and g are, and (/ o g)' = g'f o g is strictly increasing (decreasing) if 
both f',g' are strictly increasing (decreasing) and g is strictly increasing. 

• fg is strictly increasing (decreasing), hence a bijection, if /, g are both strictly increasing (decreasing), while 
(/<?)' = fg' + f'g is strictly increasing (decreasing) if f,g,f',g' are all strictly increasing (decreasing). 

□ 



5 



We have now offered a quite extensive description of the class of piecewise continuously differentiable bijec- 
tions on [0, 1] with the Costas property, and an exact characterization of the continuously differentiable bijections 
with the Costas property. What about discontinuous bijections, though? By interpreting discontinuity in the 
most extreme way, we are led back to the idea of fractals. 



5 Costas fractals 

In what follows, we establish a connection between discrete and continuum Costas permutations: we use discrete 
Costas permutations to build continuum ones through a process of multiscale rearrangement of subintervals of 
[0,1]; in other words, we build a "Costas fractal". At this moment, however, we are unable to prove the 
correctness of our construction below under the usual laws of arithmetic: we will need the equivalent of "xor" 
addition (and subtraction), namely addition without carry, in representations over an arbitrary basis. 
We will need first of all the slightly stronger definition given below: 

Definition 6. Consider a bijection [n]; / is a modulo Costas permutation iff the multiset {(i — j, f(i) — 

f(j) mod (n + 1)) : < j < i < n} is actually a set, namely all of its elements are distinct. 

Remark 4. Note that both the Golomb and the Welch constructions actually lead to modulo Costas permu- 
tations [HE]. 

oo oo 

Definition 7. Let the numbers x,y £ [0,1] be expanded over basis n £ N*: x = ~^^XiU~' , y = Vin~ l , 

i = l i = l 

where Vi £ N*, Xi, yi £ [n]. Then, we define the "no carry" addition and subtraction as: 

(xi + yi) mod n (xt - yi) mod n 

*®y = l^ - t . ^e^L ni 

Theorem 6. Let ngN and let /; : [n] — > [n], i £ N* be a sequence of (not necessarily distinct) modulo Costas 
permutations. Define a function F : [0, 1] — > [0, 1] by the following formula: 



/ oo \ oo 

\i=l / i = l 



where Vi £ N*, a; £ [n], and so that there exists no N £ N* : ai = n — 1 for i > N , unless N — 1. Then, F has 
the Costas property, when subtraction is interpreted as in Definition^ 

Remark 5. The explicit exclusion of sequences {ai}°^ 1 so that 3N £ N* : at = n — 1 for i > N is necessary in 
order to ensure that every number in [0, 1) can be expressed over base n in a unique way, otherwise some numbers 
can have 2 different expansions: a familiar example over base 10 would be that 0.5 = 0.5000 . . . = 0.4999 .... 

°° — 1 

However, we still need to represent 1 = >^ : — , hence the exception for N = 1. 

i=i 

Proof. Select 4 points in [0, 1], say x, y, x + d and y + d so that y < x and d > 0; notice that these can actually 
be 3 equidistant points if y + d = x. We need to test whether F(x) F(x + d) = F(y) F(y + d) necessarily 
implies d — 0. 



{% i + 1 
h-i = — , 
n n 

so that Vz £ [n], F(Ii ; i) = I^juy We distinguish the following cases: 



iS [n-l]}U 



il 



S~L 

1. y + d ^ x and the 4 chosen points all lie in different subintervals: then, we can write F(x) = \- ei, 

n 

^(j/) = — + £ 2, Ffx + d) = — + e 3 , and F(y + d) = — + e 4 , with Sj £ [n], 6* < -, i = 1,2,3,4. It follows 
n n n n 

\ ,% (si — S3) mod n , , (S2 — S4) mod n . . 

that F(x) F(x + d) = ^ ^ + (ei e 3 ), and F (y 0F y + d = ^ + e 2 £4 , 

n n 
where, if we assume d > 0, (si — S3) mod n 7^ (S2 — S4) mod n, by the modulo Costas property of /1, 

while |(ei £3) (e 2 £4)! < — ■ Hence, F(x) F(x + d) / F(y) F(y + d) and the proof is complete for 

n 

this case. 

2. y + d = 2; and the 3 chosen points all lie in different subintervals: then we can repeat verbatim the previous 
argument with 3 instead of 4 points. 
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3. y + d 7^ x and one pair of the 4 chosen points lie in the same subinterval, while the remaining pair lie in 
different subintervals: then, without loss of generality, assume that x and x + d lie in the same subinterval. 
In terms of the previous argument, (si — S3) mod n = 7^ (sa — S4) mod n and the proof follows again. 

4. y + d = x and the 3 chosen points lie in 2 different subintervals: then, exactly 2 points lie in the same 
subinterval, and, without loss of generality, assume they are y and y + d — x. In terms of the previous 
argument, S4 = Si, (si — S3) mod n 7^ (S2 — si) mod n = and the proof follows again. 

5. Either y + d 7^ a; and the 4 chosen points lie pairwise in the same subintervals, or y + d = x and the 3 
chosen points all lie in the same subinterval: then, assume, without loss of generality, that x and x + d lie 
in the same subinterval, and so do y and y + d. It follows that (si — S3) mod n — = (S2 — S4) mod n 
and the argument fails. 

In the last case where the argument fails, we need to refine our subinterval division. We already saw the 
first level of this division. At level k £ N, we consider the collection of intervals 



Ik;i 1 ,...,i k 



j=i 3=1 



n — 1 , n — 1 



nr 



With respect to the newly defined levels of subintervals, there are 2 possibilities: 

• The chosen points fall in a case other than for the first time in level k: then, it must be the case that: 



{fj{x3+dj)-fj(xj)) mod n _ (f k (x k + d k ) - fkjxkj) mod n 

E fj ( (Vj +dj)- fj (yj ) ) mod n = (fk(yk + d k )-f k (y k )) mod n 
ni n k 

3=1 

due to the modulo Costas property of whence F(x) Q F(x + d) 7^ ^(j/) f (2/ + d) for d > 0. 

• Otherwise, we need to consider the levels beyond level k. 

But the length of the subintervals in level k is n~ k which decays to as k — > 00; therefore, any specific selection 
of points can remain in case [5] for a finite number of levels only. This completes the proof. □ 

It is easy to see where our proof fails under ordinary arithmetic: revisiting case[T] we would need to show that, 
under the assumption that si — S3 7^ S2 — S4, which holds because /1 is a Costas permutation (we no longer need 

it to be a modulo Costas permutation), h(ei — £3) (£2 — 64) holds. Since ej < — , i = 1, 2, 3, 4, 

n n n 

1 2 
it follows that lei — £3!, lea — £4] < — and I (ei — £3) — [e% — £4)] < — , so that, if |(si — S3) — (S2 — S4)| = 1, it 

n n 
Si — S3 S2 — S4 

may still be the case that h (fi — £3) = h (£2 — £4) F(x) — F(x + d) = F(y) — F(y + d) when 

n n 
d > 0, and the Costas property fails. 

The key feature of the arithmetic proposed in Definition [7] that allowed the proof of Theorem |H] to complete 

successfully was that if, at any level of interval subdivision, the 4 chosen points were found to lie into distinct 

subintervals, the defining inequality of the Costas property would be satisfied for the chosen points. There are 

alternative arithmetics with this property: 

00 00 

Definition 8. Let the numbers x,y £ [0,1] be expanded over basis n £ N*: x = ^^a^n -1 , y = yjn~* , 

where Vi £ N*,Xi,yi £ [n]. Then, we define the "contracted" subtraction as: 

00 

. Xi - y l 

x e y = > — 5— r 

t=l 

Theorem 7. Let n £ N and let fi : [n] — * [n], i £ N* 6e a sequence of (not necessarily distinct) Costas 
permutations. Define a function F : [0, 1] — ► [0, 1] by the following formula: 



/ 00 \ 00 

F I ^2 = fi(ai)n~ l 

\i-l J i=l 
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where Vi £ N*, a* £ [n], ami so £/iai there exists no N E N* : ai = n — 1 /or i > N , unless N — 1. Then, F has 
the Costas property, when subtraction is interpreted as in Definition^ 

Proof. This is a verbatim repetition of the proof of Theorem [|J] □ 

Is it likely that Theorem [|J] still hold true for ordinary arithmetic despite the fact that our proof does not 
carry through? At this time we have no reason to believe that it does. It may still be possible to use discrete 
Costas permutations to generate a Costas fractal in the continuum, but the actual mechanism should most 
probably be different. 

6 A limiting process 

Assuming Artin's Conjecture holds true [9], which would be the case if the Generalized Riemann Hypothesis 
holds true, for any non-square integer k £ N*, k > 1 there exists an infinite sequence of primes, say {p„}„gN*, 
for which k is a primitive root. We can construct then the sequence of Welch Costas permutations corresponding 
to the primes of the sequence and the primitive root k: 

/„ : [p« - 1] + 1 -> [p n - 1] + 1, fn(i) = k 1 - 1 mod p n , i E [p„ - 1] + 1, n £ N* 

The key observation is that Vm £ N*,3iV £ N* : Vi = 1, . . . , m, f n (i) = k 1 ' 1 ; in particular, N is the smallest 
integer for which pn > k m ~ l . In other words, for any fixed number of terms, all functions of the sequence, after 
skipping a finite number of functions, have these initial terms in common. Define then the pointwise intermediate 
limit of {/njngN* to be as follows: for a fixed i £ N*, 

f(i) — lim/ n (i) = lim f n (i) := limfc 1-1 , where N is the smallest integer such that pjv > 

n> N 

in — 1 

Choose now a sequence {i„}„gN* of integers such that lim = x. We define the limit of {/nineN* 

Pn 

evaluated on {in}n€N* to be a continuum function on [0, 1] as follows: 

s(x) = lim (/(;„)) 3^ = k x , as € [0,1) 

k x — 1 

We can bring the range of s within [0, 1) as well after a linear transformation, and create: S(x) — — — ; this 

k — 1 

is the second example function in Corollary [2] 

To sum up, in the special case of an infinite sequence of Welch Costas permutations generated by a common 
primitive root k, we were able to carry out a limiting process and construct a continuum Costas permutation, 
using the property that all the members of this sequence (except possibly some of the first ones) have a common 
beginning. The limit we obtained, however, is a smooth function and not a fractal, as one might expect given 
the way Welch Costas permutations look like. 

7 Costas bijections in the rational continuum 

The idea of fractals with the Costas property in the (real) continuum was explored above in Section O where 
we saw that their implementation required special considerations. We return to this issue here, but this time in 
the context of the rationals Q = Q PI [0, 1]: in many ways the rationals stand midway between the integers and 
the reals, in the sense that they form a dense set (like the reals), but still enumerable (like the integers). We 
are about to see that these 2 properties allow us to make further progress in the subject. 

Note that Costas permutations on the rational continuum is a genuinely new problem, and in no way a 
special case of the constructions in the real continuum; the reason is that the constructions of Section U] do 
not map bijectively the rationals onto the rationals. For example, f(x) = x 2 is not a bijection over Q, as 

$x e Q : f(x) = -, say. 

The relevant definitions of the Costas property on rational bijections closely parallel the ones in Section ^. II 
(regarding the real continuum) and will not be repeated here. 
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7.1 An existence result 



In this section we offer an algorithm of considerable generality for the construction of bijections on Q with the 
Costas property. Let us begin by reordering the elements of Q as follows: we order firstly by the magnitude 
of the denominator, and secondly by the magnitude of the numerator (both in an increasing way). Explicitly, 

first come those rational numbers in [0, 1] whose denominator is 1, namely = — and 1 = —; then, those whose 

1 12 
denominator is 2, namely — ; then, those whose denominator is 3, namely — and — etc. Hence, the sequence 

looks like this: 

11213123415 
' '2'3'3'4'4'5'5'5'5'6'6"" 
Notice that the numerators are always taken to be relatively prime to the denominators in order to avoid 
duplicate entries. We denote Q equipped with this particular ordering by Qx, and its elements, in the order 
dictated by the ordering, by aio, x\, Xz, ■ ■ ■■ This ordering has the advantage that each rational is preceded 
by a finite number of rationals only (in set theoretic terminology, it does not contain any transfinite points). 
Similarly, we denote by Qy the set Q equipped with any arbitrary but fixed ordering without transfinite points, 
and we denote its elements, in the order dictated by its ordering, by yo, y\, 3/2, ■ ■ ■■ 

Consider now the following algorithm for the construction of a mapping / : Q — > Q: 

Algorithm 1. 

Initialization Choose f(xo) = yo; set Q' Y <— Qy — {yo}, Q'x <~ Qx — {%o}, X <— {a;o}, Y *— {yo}, and 
D^{}. 

Find x for y: Set Qx,nv <— Q'x, x ^ inf Qx.av, y *— mf Q' Y ; while the (multi)set {sgn(a; / — x)(x' — x, f(x') — y) : 
x' € X} U D is actually a multiset, set Qx,a,v <— Qx,a,v — {%}, x <— inf Qx,&v, and repeat. Set f(x) — y, 
D «- {sgn(x' - x)(x' - x, f(x') - y) : x' € X} U D, Q' Y «- Q' y - {y}, Q' x <- Q' x - {x}, X <- X U {x}, 
Y<-Yu{y}. 

Find y for x: Set Qy,hv <— Q'y , V ^~ m ^ Qy.hv, x <— inf Q' x ; while the (multi)set {sgn(a;' — a;) (a;' — x, f(x') — y) : 
x' € X} U D is actually a multiset, set Qy, av <— Qy, a v — {y}, y *— inf Qy. av , and repeat. Set f(x) = y, 
D «- {sgn(x' - x)(x' - x,f(x') - y) : x' G X} U D, Q' Y «- Q'y - {y}, Q'x <- Q'x - {*}, X <- IU {x}, 
Y^YU{y}. 

The algorithm needs to be supplied with a step sequence before execution begins. For the purposes of the 
correctness proof the exact step sequence is unimportant (this is yet another degree of freedom of the algorithm), 
as long as the following rules are observed: 

• Initialization is run first and only once; 

• Neither Find x for y nor Find y for x is run infinitely many times in a row. 



For example, when Qy = Qx and the steps are run alternatingly, we get /(0) = 0, /(l) = 1, / 

I- f {k) = * / (i) = f etc - 

Theorem 8. Algorithm]!] produces infinitely many bijections f : Q — > Q with the Costas property. 

Proof. In order to prove the correctness of Algorithm [T] above, we need to demonstrate that a) Vy £ Qy, 3!x € 
Qx : f(x) = y, and b) Va; £ Qx,3y £ Qy '■ f(x) = y. To begin with, note that the construction algorithm 
above guarantees that the constructed / has the Costas property and that every y £ Qy appears in the range of 
/ at most once. We only need to show that the algorithm never gets "stuck" , namely that the two while loops 
always exit. 

• For a given x, is it possible to assign a value to f(x)? In other words, if A C Q'y is the set of all values 
f(x) can take without violating the Costas property of /, is it true that A 7^ 0? The answer is in the 
affirmative, as, intuitively, we can see that the Costas property restrictions impose only a finite number of 
constraints on f(xi), while Q' Y is countably infinite. Rigorously, we have to check 2 conditions: 

— Let A\ C Q'y be the set of possible values for f(x) for which sgn(a;' — x)(x' — x, f(x') — f(x)) = sgn(a;' — 

x")(x' — x",f(x') — fix")) is never true for x'.x" G X. We show that A\ 7^ 0. In fact, consider -, 

V 

where p is a prime that does not appear as a factor in the denominator of some f(x'), x' G X: choosing 
f(x) = — , it follows that f( x ') contains p as a factor in the denominator, while f(x') — f(x") 
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does not, hence they cannot be equal, and therefore that — G A\ 7^ as promised. Clearly, there are 

P 

infinitely many choices for p possible, so Ai contains actually infinitely many elements. 

— Let A2 C Ai be the set of possible values for fix) for which sgn(a;' — x)(x' — x,f(x') — f(x)) = 
sgn(x" — x)(x" — x, f(x") — f(x)) is never true for x' , x" G X. We show that A2 7^ 0. In order for one 
of these equalities to hold, x must be the midpoint of x' and x" , while at the same time f(x) be the 

midpoint of fix') and f(x"). Choosing fix) — — where p is as above, and writing x = — , x" = — , 

p Vi V2 

we need to investigate whether the following is possible: 

+ — J = (ui.vi) = {U2,V 2 ) = 1, p Jfvi,v 2 . 
2 \vi v 2 ) p 

This implies p(«i«2 + U2V1) = 2v\V2, and therefore that p\2v\V2 =>■ p\2 p — 2. Hence, A2 does 

contain all points of the form -, too, where p does not divide the denominator of some fix'), x G X 

P 

(which are infinitely many), except possibly for — ; in any case A2 7^ 0. 

But A2 = A, hence A 7^ 0, a contradiction; therefore, f(x) can assume a value without / losing the Costas 
property. 

• For a given y, is it possible to find x G Q'x so that fix) = yl In other words, if A C Q'x is the set of all 
values x for which fix) can be y without violating the Costas property of /, is it true that A 7^ 0? The 
answer is in the affirmative as well, and the argument is an almost verbatim repetition of the argument 
above. Rigorously, we have to check 2 conditions: 

— Let A\ C Q'x be the set of possible values for x for which sgn(a;' — x)ix' — x, fix') — y) = sgn(a; / — 

x")(x' — x",f(x') — fix")) is never true for x',x" G X. We show that Ai 7^ 0. In fact, consider -, 

P 

where p is a prime that does not appear as a factor in the denominator of some x' G X: choosing 

x = -, it follows that x contains p as a factor in the denominator, while x' — x" does not, hence 

P P 

they cannot be equal, and therefore that - G Ai 7^ as promised. Clearly, there are infinitely many 

P 

choices for p possible, so A\ contains actually infinitely many elements. 

— Let A% C A\ be the set of possible values for x for which sgn(a;' — x)(x' — x, f(x') — y) = sgn(a;" — 
x)ix" — x,fix") — y) is never true for x',x" G X. We show that A2 7^ 0. In order for one of these 
equalities to hold, x must be the midpoint of x' and x" , while at the same time y be the midpoint of 

fix') and fix"). Choosing x = — where p is as above, and writing x' = — , x" = — , we need to 

p Vl V2 

investigate whether the following is possible: 

( — + — j = I = (u 2 ,V 2 ) = 1, P J(Vl,V2. 

2 \vi v 2 J p 

This implies p(uiU2 + U2V1) — 2viv 2 , and therefore that p\2viV2 p\2 p = 2. Hence, A2 does 
contain all points of the form -, too, where p does not divide the denominator of some fix'), x G X 

(which are infinitely many), except possibly for i; in any case A2 7^ 0. 

But A2 = A, hence A 7^ 0, a contradiction; therefore, there exists a x : /(a;) = y without / losing the 
Costas property. 

This completes the proof. □ 

Remark 6. Intuitively, the mechanism responsible for the flexibility of the algorithm is the opportunity the 
countable infinity of the rationals offers for "double deference of all difficulties for a future time": when faced 
with the difficulty of assigning a value to / at a given point, we always have infinitely many possibilities, out of 
which some will work; this in turn creates the difficulty of assigning the values we skipped to some point, but, 
when faced with this difficulty, we again have infinitely many points waiting for an assignment, out of which 
some again will work; but in choosing one we once more skip some points, and we need to choose values for 
them, hence the cycle restarts. 

This interplay is precisely what we cannot do with a finite set, hence the contrast between the easiness of the 
Costas construction over the rationals, as opposed to the intractability of the classical construction of Costas 
arrays. 
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Remark 7. The above proof makes heavy use of the enumerability of the rationals, and therefore cannot be 
readily extended to the reals, who lack this property. 

It may come as a surprise that we can extend the algorithm even further: 

Theorem 9. Algorithm\7\ will produce a bijection f : Q — > Q with the Costas property even if one of the steps 
From x to y or From y to x is applied infinitely many times in a row. 

Proof. Let us consider the case where From x to y is run infinitely many times in a row immediately after 
Initialization. This causes no loss of generality: the case where From y to x is run infinitely many times in a row 
immediately after Initialization is completely dual (observe the duality in the proof of Theorem [8| , while the 
more general situation where finitely many alternations between the 2 steps occur before the algorithm "locks" 
in one can be considered to fall within one of the 2 cases we just mentioned, but with a different, more extensive 
Initialization. 

Assume then that we go through x £ Qx one after another and we try to assign values to f{x) £ Qy while 
retaining the Costas property. The proof of Theorem [H] guarantees that we will succeed for all points. What 
we need to worry about is whether some y £ Qy will be left out in the process: in other words, we know that 
Vx £ Qx, 3y £ Qy : f(x) = y, but we still need to know that Vj/ £ Qy, 3!x £ Qx : f(x) = y. 

Assume then that at some step of the algorithm we find that y £ Q' Y has been skipped, and is the smallest 
element of Qy that has been skipped. Will the algorithm ever "pick it up"? As before, let us denote by A C Q'x 
the set of all available x for which we can set f(x) — y without violating the Costas property; we need to show 
that A ^ |. Because we proceed through Qx sequentially from the beginning, at the particular step of the 
algorithm we find ourselves there exists xo £ Qx '■ X — {x £ Qx ■ x < xo} (remember that < refers to the 
ordering of Qx, not the usual ordering!). 

Consider aid Q' x which is of the form -, p prime, say x- As in the proof of Theorem [5] we need to show 

P 

2 things: 

• sgn(a;' — x)( x ' ~ X,f( x ') ~ V) = sgn(a;' — x")(x' — x" , f(x') — f(x")) is never true for x',x" < \. The 
additional complication here is that at the current step of the algorithm we know the values of / up to xo, 
but we endeavor to prove a property that holds for x < x, i- e - involving future values! The way to avoid 
the complication is to apply our favorite argument on the first coordinate only, disregarding entirely what 
the values of / are: x — x" cannot contain p as a factor in its denominator, while x~ x ' does, hence they 
cannot be equal. It follows that x wm belong in A as long as it satisfies the second condition we are now 
about to test, and also that x can actually be chosen among infinitely many points. 

• sgn(x' — x)( x ' ~ X,f{ x ') ~ u) = sgn(x" — x)( x " ~ X,.f( x ") ~ v) ls never true for x',x" < x- In order to 
check this we repeat verbatim the proof of Theorem [5] we assume that x ls the midpoint of some x' and 

x" , and then show this is impossible, unless perhaps x= ~- It follows that - satisfies this condition too, 

2 p 

with the possible exception of when p — 2. But this still leaves infinitely many points of the form -, p 

P 

prime, in A, hence in particular A 7^ 0. 
This completes the proof. □ 

7.2 An explicit construction 

Algorithm [T] is not exactly constructive; we cannot, for example, readily compute what / ( iq2§ ^ 1S e< l ua ^ 

We propose here a constructive algorithm for the construction of a Costas permutation on the rationals; the 
catch is, however, that it only works on a subset of Q. 

Definition 9. We define the set of prime rationals Qp in [0, 1] to be the subset of Q with prime denominators; 

namely Qp = < — : i € [p — 1], p prime > . 

IP J 

Theorem 10. For each prime p, consider a Welch Costas permutation f v : [p — 1] + 1 — > [p — 1] + 1 constructed 
in ¥(p), and consider the set of points S(p) = < (—, j : i £ [p — 1] + ll. The set S = S(p) is a 

p prime 

Costas permutation on Qp. 

Proof. S is clearly a permutation. We need to show that the distance vectors between all pairs of points are 
distinct. 
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• Choose 4 points in the same S(p): the Costas property of f p guarantees the 2 distance vectors they define 
are distinct. 

• Choose 2 points in S(p) and 2 points in S(q), q 7^ p: the first distance vector has coordinates that are 
fractions over p, while the second over q, hence they cannot be equal. 

• Choose 2 points in S(p), a point in S(q), and a point in S(r), where p,q,r are distinct primes: the first 
distance vector has coordinates that are fractions over p, while the second over qr, hence they cannot be 
equal. 

• Choose a point in S(p), a point in S(q), a point in S(r), and a point in S(s), where p,q,r,s are distinct 
primes: the first distance vector has coordinates that are fractions over pq, while the second over rs, hence 
they cannot be equal. 

This completes the proof. □ 

8 Conclusion 

In this work, we have made 4 main and original contributions to the subject of Costas arrays: 

• We defined the Costas property on a real continuum function in 2 ways, through distance vectors between 
points and through the autocorrelation, and we showed that the 2 definitions are equivalent. We also 
showed that real continuum Costas bijections can be used in the same applications as discrete Costas 
arrays, by designing signals with the appropriate instantaneous frequency, which has been made possible 
by the recent advances in the field. Subsequently, we studied similarly the Costas property on rational 
continuum functions. Essentially, we have now translated the entire framework of Costas arrays in the 
continuum. 

• We showed that real continuum Costas bijections exist and we offered some examples; we characterized 
completely the continuously differentiable Costas bijections in terms of the monotonicity of their derivative, 
and we also obtained some good results for the case where the bijections are only piecewise continuously 
differentiable. 

• We investigated whether it is possible to construct fractal bijections with the Costas property, perhaps by 
employing discrete Costas arrays as building blocks. We answered that in the affirmative under nonstan- 
dard arithmetic laws (where addition and subtraction take place without carry, or where the contribution 
of the least significant digits of the points to their distance is deemphasized) in the real continuum; under 
ordinary arithmetic we have no reason to believe that the result still holds true. 

• We proposed a very general and flexible algorithm for the construction of Costas permutations over the 
rationals, that is not, however, entirely constructive. We were also able to formulate such a constructive 
algorithm, but its applicability is limited over a subset of the rationals. 

Overall, it comes to us as a surprise that it was relatively simple to construct smooth continuum functions 
with the Costas property, whereas all efforts to create a fractal Costas real bijection were unsuccessful (under 
ordinary arithmetic). Intuitively, given the irregularity of discrete Costas arrays, we would expect the known 
construction methods for Costas arrays to generalize in a natural way in the real continuum leading to a fractal; 
however, a direct recursion, such as our attempt in Section seems to be inappropriate, unless we change the 
arithmetic we use. It may still be possible to construct a Costas fractal bijection based on discrete Costas arrays 
through a different, less obvious mechanism, and we challenge the reader to discover such a mechanism. 
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